Averages and critical exponents in type-Ill intermittent chaos 
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The natural measure in a map with type-Ill intermittent chaos is used to define critical exponents 
for the average of a variable from a dynamical system near bifurcation. Numerical experiments were 
done with maps and verify the analytical predictions. Physical experiments to test the usefulness 
of such exponents to characterize the nonlinearity at bifurcations were done in a driven electronic 
circuit with diode as nonlinear element. Two critical exponents were measured: v = 0.55 for the 
critical exponent of the average of the voltage across the diode and j3 = 0.62 for the exponent of 
the average length of the laminar phases. Both values are consistent with the predictions of a type- 
Ill intermittency of cubic nonlinearity. The averages of variables in intermittent chaotic systems 
is a technique complementary to the measurements of laminar phase histograms, to identify the 
nonlinear mechanisms. The averages exponents may have a broad application in ultrafast chaotic 
phenomena. 
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I. INTRODUCTION 

Intermittent chaos is the phenomenon shown by sys- 
tems exhibiting long sequences of periodic-like behaviour, 
the laminar phases, separated by comparatively short 
chaotic eruptions. Intermittent chaotic systems have 
been extensively studied since the original proposals by 
Pomeau and Manneville 0, H H || classifying type-I, 
II, and III instabilities when the Floquet multipliers of 
a map crosses the unity circle. Although other mech- 
anisms may occur leading to intermittency these three 
cases are the most simple and most frequently encoun- 
tered in low dimensional systems. The class of intermit- 
tent chaos studied by Pomeau and Manneville includes 
the tangent bifurcation, leading to intermittency of type 
I, when the Floquet 's Multiplier for the associated map 
crosses the circle of complex numbers with unitary norm 
through +1; the Hopf bifurcation, leading to type-II in- 
termittency, which appears as two complex eigenvalues 
of the Floquet 's Matrix cross the unitary circle off the 
real axis; and the sub-critical period doubling, leading 
to type-Ill intermittency, whose critical Floquet 's Multi- 
plier is —1. 

Many experimental evidences for these intermittent 
behavior have appeared in the literature. The type- 
Ill has been reported for electronic nonlinear devices 
H10, lasers [g , and biological tissues . A signature 
for this intermittency is given by the critical exponent 
describing the dependence of average length of nearly 
periodic phases, that is, laminar phases, with the con- 
trol parameter. Histograms of number of laminar phases 
longer than a given duration are related to the normal 
form nonlinearity describing the system 0, EU ■ Different 
nonlinear power in the model imply different exponents, 



as described by Kodama et al. 01 • Kim et al. 3 studied 
the exponent for the average length of laminar phases as 
function of the reinjection probability. Herein another 
average is explored to characterize the intermittency and 
its nonlinearity. It is the average of one variable of the 
system. Near bifurcation, approximate expressions for 
the natural measure, or probability density, are obtained 
for a map with type-Ill intermittency, in section^ Criti- 
cal exponents are defined from analytical approximations 
and shown to be useful in characterizing the bifurcation. 

Numerical experiments with the maps are presented in 
section ITTT1 verifying the proposed exponents. To test in 
a real physical system, a nonlinear circuit with a diode, 
similar to the one used in the early demonstrations of 
chaotic universal properties 0, 0, E2> E3 j was set up as 
described in section Hvl and used to verify the exponents. 
The voltage across the diode, which is a dynamical vari- 
able of the system, was simultaneously measured in time 
series and average. The type-Ill intermittency bifurca- 
tion is well characterized both, using the experimental 
next peak value map and the average. 



II. NORMAL FORM MAP WITH TYPE-III 
INTERMITTENCY 

To establish the new critical exponents for the averages 
one begins with the normal form of the map that has 
type-Ill intermittency 



M(X) = -(1 + e')X + aX 2 + rj'X 3 



(1) 



The bifurcation, when X = ceases to be a stable fixed 
point, occurs at e' = 0. The second application of this 
map, in the approximation of small e' and X is given by 



M 2 (X) = (1 + 2e')X + i]X ? ' 



(2) 



where rj 



-2(7/ 



If this coefficient is posi- 



tive and e' > one has type-Ill intermittency around 



2 



X = 0, provided a reinjection mechanism is introduced 
in eqs. Q and (01. When a^O the map of eq. £[J is 
nonsymmetrical in X. Generally, any odd exponent in 
the nonlinear term of eq. J2J leads to the intermittency. 
Thus, a map (second iterate) defined in the interval [0, 1] 
with the type-Ill intermittency is 3] 



X n+1 = [(1 + e)X n + X*] mod 1 



(3) 



The value of z > 3 describes the nonlinear dependence 
of the subcritical bifurcation at e = 0. For e <C 1 the 
fixed point X = is unstable, but many iterates of the 
map, shown in Figure^ falls near zero. These are called 
the laminar phase iterates. 
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FIG. 1: Map with the normal form that gives type-Ill inter- 
mittency. The value z — 3 was used in eq. ®- 

For small e the natural measure is obtained, following 
the steps of Manneville [14J , as 



fi(e, z, X, dX) 



dX 



dX 



e + X^-V 



(4) 



Using this measure, the average of X is 



(X(e,z)) = 



dX 



+ X(*-V 



X dX 

+ x( z -v 



(5) 



Its dependence on e and z can be analytically established 
for small e. 

The simple case z = 3 reduces to 



(X(e,z = 3)) «-i e 1 /2 lne 



For z > 3 the general asymptotic expression is 
(X(e,z)) oc e u 
where v is the new exponent, whose value is 

1/ = 1/(2-1) 



(6) 



(7) 



(8) 



This is the exponent that can be extracted from numer- 
ical and experimental systems to obtain the value of z. 
Similar expressions can be obtained for the second iter- 
ated map with negative values of X. If the original first 
iterated map, giving eq. © is fully odd in X the total 
average must be zero. Once the reinjection or the map is 
not symmetrical, i. e., the value of some even nonlinear 
power (less than z) is non zero, the exponents of the av- 
erages can be obtained from the second iterate maps as 
eq. Q. A detailed study of these symmetry properties 
will be published elsewhere [15| . One must recall that, 
in type-Ill intermittent chaos, an exponent already ex- 
ists in the literature for the average length of laminar 
phases, given by 013 P = (z- 2)/{z ~ 1). The simple 
relation exists, giving v + /3 = 1. Therefore, as z varies 
through 3, 5, . . . the relative variation of these exponents 
are /S.v/v > A(3/[3. This means equal or better sensi- 
tivity in the determination of v with respect to (3. The 
most important fact is that those two exponents should 
be searched for independently in any experiment. 



III. NUMERICAL AVERAGES AND 
EXPONENTS IN THE MAP 

To test the predictions of eq. JBJ numerically the map 
of eq. 0) was restricted to \X\ < 1 by a modulo one 
reinjection: When the iterate gives 1 > 1 or 1 < - 1 
the integer part is subtracted. Averages were obtained 
directly from iterates of the map starting from random 
initial condition. A transient of 2 x 10 iterates was elim- 
inated from a total of 10 6 calculated values at each of 10 3 
steps of e', between zero and 1CP 3 . For the average to be 
non null the value of a has to be different of zero. This 
breaks the symmetry between positive and negative val- 
ues of X. Once the average is non null, its behavior for 
small e follows eq. ©■ This is shown in Figure [3 The 
a = case, where the average is always zero, is not repre- 
sented. The same slope is obtained for the averages with 
different values of a and Tj'. The resulting averages have 
the same dependence in e as the one obtained from the 
map of second iterates (eq. ©), calculated with the same 
number of iterates; all coinciding with the expression of 
eq. ©. 

Figure [3] shows numerically calculated averages ob- 
tained directly from iterates of the map eq. JSJ, when 
z = 3 and z = 7. Again the initial condition at each 
value of e was taken at random. A transient of 10 5 it- 
erates was eliminated from a total of 2 x 10 7 calculated 
values at each of 10 3 steps of e, between zero and lfU 1 . 

For e < 10~ 2 the z = 3 and 7 results are very 
well superimposed by the analytical curves made from 
equation 10, with v = 0.50, and equation J7J, with 
v = 1/6 « 0.17, respectively. The two behaviors are 
clearly distinguished in the log-log plot. For the purpose 
of comparison with the already established exponents one 
should notice that the exponents for the average laminar 
phase are given in [lj} as (z — 2)/(z — 1) that is (3 — 0.50 
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FIG. 2: Average of X from the map in equation Q, with the 
modulo 1 reinjection and calculated with 10 6 iterates at each 
of 10 3 steps of the control parameter e. The log-log plot shows 
the same slope for all averages, independent of the value of a 
and r\ . 
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FIG. 3: Average of X from the map in equation Q, calcu- 
lated with 2 x 10 7 iterates at each of 10 3 steps of the control 
parameter e. White squares are used for z = 7 and filled 
circles for z = 3. The thin lines correspond to the predicted 
values given by eqs. © and 10. 



and (3 = 0.83, respectively, verifying v + (3 = 1. 

Up to this point all results have concerned discrete 
maps. Dynamical fluxes with associated maps having in- 
termittent chaos show similar critical exponents for the 
average of its continuous variables. This has been ver- 
ified numerically with tangent bifurcation of type-I in- 
termittency [TtI Il8j and will be demonstrated here for a 
physical chaotic oscillator. For the type-Ill intermittency 
the map extracted from a flux has the form of eq. J3J for 
its second iterates. A discretely sampled continuous flux 
may be viewed as a finite set of discrete stroboscopic 
maps, all with the same stroboscopic frequency but stro- 
boscopic phases ranging from to 2tt in uniformly sepa- 
rated steps. The number of such maps gives the ratio of 
sampling frequency to stroboscopic frequency. The time 
average of the continuous flux is the average of the time 
averages of these maps. If a map from this series is non- 
symmetrical, only in very special cases it would happen 



that this asymmetry be exactly canceled by the summa- 
tion on the other maps. Thus, in general, its expected 
that if the flux leads to a nonsymmetrical map (no odd 
symmetry with respect to the unstable fixed point) the 
average of its continuous variable will exhibit the same 
exponent v of eq. JSJ. 



IV. EXPERIMENTAL EXPONENTS IN A 
NONLINEAR CIRCUIT WITH 
INTERMITTENCY 



The physical experimental system to test the above 
exponents consisted of an RLC series circuit, where a 
nonlinear capacitance was implemented by a p-n junc- 
tion diode 0, 0. The circuit, presented in figure 0] 
was driven by a frequency and amplitude controllable 
oscillator with output impedance of 50 Ohms. An induc- 
tor of 0.1 Henry, an external resistor of 13 Ohms and a 
(typical 1N4007) diode formed the circuit, whose linear 
oscillation regime had resonance frequency at 150 kHz. 
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FIG. 4: Experimental diagram of the RLD circuit showing 
type-Ill intermittent chaos. The voltage across the diode was 
the dynamical variable measured both in its peak value, to 
make maps, and integrated to give the averages. The circuit is 
drawn with thick lines, while thin lines are the high impedance 
probe and integrator. 

Dynamical bifurcations were produced by scanning the 
external drive frequency. Time series of the value of the 
voltage across the diode and the current in the circuit 
were collect with a 12 bits resolution converter. The sam- 
pling rate was 10 7 sample/s. Thus 200 points were saved 
on each oscillator cycle. Capturing 10 4 points at each 
value of the control frequency, a simple software searched 
for the maxima in the series. In wide range scans the typ- 
ical bifurcation diagrams given by the peak value of the 
voltage across the diode, show clearly the well known re- 
sults of period doubling cascades, chaotic windows and 
tangent bifurcations from chaos into periodic windows. 

lamia. 

A specific bifurcation, with intermittency and no bista- 
bility, was found and studied, scanning the external oscil- 
lator frequency from 48708 Hz to 48768 Hz by 300 equal 
small steps. The drive voltage amplitude was fixed at 
2.7 V. The circuit was mounted inside an isolating box 
to prevent thermal drift effects 0|. The value of a con- 
trol parameter e is obtained as the frequency detuning 
step divided by the critical frequency of the bifurcation, 
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48762 Hz. Thus e varied in the range of 10~ 4 -10" 3 . A 
segment of the voltage pulses is shown in figure [SJ The 
signature of type-Ill intermittency is observed with the 
laminar events corresponding to the uniform oscillations 
alternating their peak value around IV. 
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FIG. 5: Voltage pulses measured across the diode in the RLD 
circuit. Successive maxima alternate their position respective 
to the 1 V value until a chaotic burst occurs and reinjects 
back to the vicinity of the unstable orbit, typical of type-Ill 
intermittency 

Multibranch maps constructed from the maxima of 
the voltage across the diode, extracted from series with 
6 x 10 5 , could be approximated to one-dimensional ones, 
expected in the limit of infinite dissipation. Second re- 
turn maps and laminar phase histograms, not shown 
here, indicate type-Ill intermittency. Figure shows the 
average length of laminar phases for different values of e. 
Each histogram used was obtained with 6 x 10 5 points. A 
theoretical fitting ^(j is best with an exponent (3 — 0.62. 
Notice that for the map the predictions are j3 = 0.5 for 
z = 3 and [3 — 0.75 for z = 5. However the experimental 
data always has excess of laminar events identified with 
short length |{J and this effect gives a bigger experimen- 
tal value for (3 in the fittings. Therefore, z = 3 is the best 
odd value. The excess of short laminar phase events may 
be related to the finite dissipation rate in the phase space 
of the experimental system and the consequent non uni- 
dimensional maps and also to a nonuniform density of 
probability for the reinjection, as discussed by Kim et 
al. ||. 

The peak value of the voltage, which gave the second 
return map and the histograms is shown in Figure Ufa) . 
The laminar phases are the oscillations with repetitive 
visits to the maximum value near IV in the figure. As 
the segments represented in figure EJa) are short in time 
for many values of parameters they show a single laminar 
event with all points accumulated around the unstable 
orbit. Also shown is the simultaneously acquired average 
of the voltage across the diode. It was obtained with a 
simple electronic integrator, having a time constant of 3 
seconds. To account for long laminar phase events and 
(which is equivalent) decrease the average fluctuations 
the scan lasted 50 minutes. 
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FIG. 6: Measured average lengths of laminar oscillations in 
the driven circuit with a diode. The thick line is the theoret- 
ical fitting for type-Ill intermittency with exponent /3 = 0.62 
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FIG. 7: (a) Peak voltage across the diode in the type-Ill bi- 
furcation of the circuit, (b) Simultaneously measured average 
voltage. The thick line is a fitting of eq. © with exponent 
v = 0.55 



The experimental average voltage was fitted to the ex- 
pression 



(X(e)) 



-e v lne 



(9) 



The exponent v = 0.55 gave an excellent agreement with 
the experimental plot, as shown in figure Attempts 
of fittings to higher values of z, i.e., to eq. Q, failed. 
It is worth noticing that the value predicted for z = 5 
is v = 0.25. Thus, the experimental average consistently 
verifies z — 3 for the nonlinearity of this bifurcation in the 
circuit. The confidence for the experimental values of f3 
and v from fittings using standard \ 2 procedure is better 
than 2%. Therefore, the result (3 + u = 1.17 ±0.02 shows 
that a discrepancy remains between the experiments and 
the unidimcnsional map model. While the bigger value 
obtained for (3 has been attributed to an excess of short 
laminar phase events [4j no such study of deviations exists 
for the exponent in the averages. 
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V. CONCLUSION 

In conclusion, critical exponents for the averages of 
one dynamical variable are established analytically for 
type-Ill intermittent chaotic maps. Those exponents are 
directly related to the nonlinear power law of the normal 
form of the maps. They have a simple relation with the 
exponents of the average length of laminar iterates in 
the same systems. All these properties were verified in 
numerical experiments with maps. 

Physical experiments with a continuous flux were also 
done to demonstrate the exponents in averages. The av- 
erage voltage across a diode in a chaotic electronic cir- 
cuit was measured while the drive frequency was scanned 
through bifurcations. A bifurcation from chaos into pe- 
riodic pulsation, shown to be type-Ill intermittent, gives 
an exponent in agreement with the cubic nonlinearity. 
This result is consistently verified in the exponents of 
the average length of laminar phases, extracted from his- 
tograms of the peak pulse voltages in the circuit. 

Averages of dynamical variables have been proposed 
to get the signature of the Lorenz chaos bifurcation 0] , 
have been numerically studied in critical bifurcations 
|20|. and experimentally measured in bifurcating pulsed 
lasers |2l|. However no systematic study of critical ex- 
ponents have been done. 

The technique of measuring averages of dynamical 
variable in intermittent chaos is a complementary pro- 
cedure to investigate bifurcations of nonlinear systems. 
The experimental average may also be advantageous 
when detection noise for an specific variable has a 
bandwidth overlapping the frequency bandwidth of the 



chaotic oscillations. For systems with high frequency 
noise it is naturally bound to be more sensitive. Its ex- 
perimental motivation is therefore enhanced to charac- 
terize ultrafast chaotic oscillators, like diode lasers, using 
slow time detection techniques. One extension relevant 
to the work presented here is the study of the expo- 
nents for the averages in bifurcations with nonuniform 
reinjection in the intermittency, as studied by Kim et 
al. |4| . A relation between exponents of the averages and 
the exponents of the average length of laminar phases 
events should exist generalizing the u + /3 = 1 given here. 
Another potential use for the averages near bifurcation 
would be to complement the confrontation between ex- 
perimental data and model as extracted from nonlinear 
data analysis [22, ■ The models inferred from the data 
analysis must have bifurcations consistent with the ex- 
perimental system. These bifurcations could be tested 
comparing both the exponents of laminar events and the 
exponents of averages of dynamical variables. The Crit- 
ical exponents have been introduced for many types of 
bifurcations in chaotic systems [24^. Their presence in 
simple, experimentally accessible, statistical properties, 
as the averages and its higher moments, are under inves- 
tigation. The earliest citation of averages in chaos can 
be traced to the original propositions of unpredictability 
in deterministic chaos [25|. 
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